function diff = delta_u_p(m,param)

    % derivative of delta(m) in the expansion region

    diff=(1-param.omega_1)./param.alpha./param.c ...
        +param.delta_2.*m.^(1/(1-param.alpha)).*(1/(1-param.alpha))./m;


end 